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ABSTRACT 

We present XMM-Newton observations of the WC binary 6 Muscae (WR 48), the second brightest Wolf-Rayet binary in optical 
wavelengths. The system consists of a short-period (19.1375 days) WC5/WC6 + 06/07V binary and possibly has an additional O 
supergiant companion (O9.5/B0Iab) which is optically identified at a separation of ~ 46 mas. Strong emission lines from highly 
ionized ions of C, O, Ne, Mg, Si, S, Ar, Ca and Fe are detected. The spectra are fitted by a multi-temperature thin-thermal plasma 
model with an interstellar absorption Nh = 2-3 xlO 21 cm 2 . Lack of nitrogen line indicates that the abundance of carbon is at least 
an order of magnitude larger than that of nitrogen. A Doppler shift of ~ 630 km s -1 is detected for the O VTJJ line, while similar 
shifts are obtained from the other lines. The reddening strongly suggests that the emission lines originated from the wind-wind shock 
zone, where the average velocity is ~ 600 km s _I . The red-shift motion is inconsistent with a scenario in which the X-rays originate 
from the wind-wind collision zone in the short-period binary, and would be evidence supporting the widely separated O supergiant 
as a companion. This may make up the collision zone be lying behind the short-period binary. In addition to the emission lines, we 
also detected the RRC (radiative recombination continuum) structure from carbon around 0.49 keV. This implies the existence of 
additional cooler plasma. 
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1. Introduction 



Wolf-Rayet (WR) stars are the evolved descendants of massive 
O stars, whose optical spectra are characterized by strong he- 
lium, carbon, nitrogen and oxygen emission lines. These arise 
from hot stellar winds with typical terminal velocities of m 1000- 
3000 km s _I and mass loss rates of the order of 10~ 5 Mq yr _1 . 
The visible spectra of WN stars and WC stars are dominated by 
nitrogen and carbon emission lines, respectively. 

Pollock ( |1987 1 used the Einstein observatory, which covered 



the 0.2-4.0 keV energy band, to discover that several WR+O 
binary systems tend to be brighter in X-rays than single WR 
stars. This was confirmed by a ROSAT survey (Pollock 1995). 
The strong hard X-ray emission from WR binaries can be inter- 
preted to arise from a collision of stellar winds from the WR star 
and that of the O star (e.g., Koyama et al. 1990[ l. Recently, |Schild| 
et al.| ( |20Q4) > discovered the RRC (radiative recombination con- 
tinuum) structure of carbon in the spectrum of the WC+O binary 
y 2 Velorum, which shows the existence of another cooler com- 
ponent. 

6 Muscae (WR 48, HD 1 13 904) is known to be a complex 
system. The history of observations and interpretations are sum- 
marized in Hill et al. (2002). Based on its spectrum, this sys- 
tem is known to be a binary or a line-of-sight double, and is 
the second brightest WR binary in optical wavelengths (Moffat 
& Seggeniss[T977| . Moffat & Seg geniss|dl977b showed that the 
WR star emission lines exhibited radial velocity variations with 
a period of 18.341 d but the absorption lines from the late O 
supergiant (O9.5/B0Iab; van der Hucht|200T i seemed to be sta- 
tionary. This leads to two possible interpretations; one is an ex- 
tremely large mass ratio, and the other is the hypothesis that the 
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Fig. 1. XMM-Newton MOS 1 image (0.3-10 keV band). 



WR component shares its 18.341-d orbit with a third hidden star, 
and then the O supergiant star is either in a much wider, slower 
orbit with the other two or is a line-of-sight coincidence. The 

~ fT987] con- 



linear polarization data obtained by |St.-Louis et al. 



firmed that the O supergiant does not participate in the short- 
period orbit. Hartk opf et al.| (fl999) resolved the O supergiant to 




Fig. 2. XMM-Newton RGS/MOS/PN spectra of 9 Muscae. The solid lines indicate the best-fit wabs * gsmooth * (redge + 2T 
vpshock) model in Table [T] 



be 46 ± 9 mas away from the short-period binary using speckle 
interferometry. 

Non-thermal emission is detected in radio wavelengths from 
this system (|Leitherer et al.l 1997 Chapman et al. 1999). 



Dougherty & Williams ( 2000| i discussed, in their study of 23 
WR stars, the possibility that the non-thermal emission arises 
not from a wind-wind collision zone in the short-period binary 
but from that between the supergiant and the short-period bi- 
nary. The discussion is based on their obtained correlation be- 
tween orbital period and radio spectral index in WR binaries. 
The data point for an orbital period of about 18 days and the 

spectral index of 9 Mus ( 0.4) is offset from the correlation, 

while if the orbital period is that estimated for the wide binary 
(at least 130 yr), the data point lies on the correlation. The cor- 
relation would be reasonable if we think that detection of radio 
non-thermal emission is difficult in short-period systems because 
the wind-collision zone is deep inside the opaque region of the 
stellar wind from the WR star. 

In X-ray wavelengths, the Einstein observatory detected 
strong X-ray emission with an obs erved lumino sity of 2.0 x 10 33 
erg s 1 in the 0.2^.0 keV band ( |Pollock|l987| , and the ROSAT 
observatory detected emission at 1.4x 10 erg s -1 in the 0.2-2.4 
keV band ( Pol lock et al. 1 1 995 1 . However, despite the past obser- 
vations, the origin of the X-ray emission is not well understood. 



The large effective area and the excellent spectral resolution over 
a broad band of XMM-Newton allows us to measure the detailed 
structure of atomic lines. We present here XMM-Newton obser- 
vations of 9 Muscae. 

We adopt the short orbital period of 19.1375 days is derived 
most rec ently (|Hill et al. 2002} and assume the distance to be 
2.27 kpc ( |van der Hucht|2001) . 



2. Observations and data reduction 

XMM-Newton was launched in 1999 December from French 
Guyana. The observatory consists of three X-ray telescopes and 
five X-ray instruments: three European Photon Imaging Cameras 
(EPIC: PN and two identical MOS CCDs), which have high sen- 
sitivity up to 10 keV, and two identical reflection grating spec- 
trometers (RGS), which have excellent energy resolution in the 
soft X-ray band (0.33-2.5 keV). 

9 Muscae was observed by XMM-Newton between 2004 July 
20 11:21:43 UT and 2004 July 21 20:36:35 UT, giving 119 692 
s of observation time. For the orbital solution of the WC star, we 
adopted that obtained in |Hill et al. (2002]!: 



HJD(af = 0) = 245 1377.51 + 19.1375£ 



(1) 
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where the eccentricity is zero and E is the number of orbits since 
the passage on HJD 245 1377.51. Phase 0.0 is defined as the 
time when the WR star occults its companion. Our observations 
cover the orbital phase from 0.596 to 0.668, during which the O 
star was located in front of the WR star. Each phase has errors 
+0.024 and +0.012. The former corresponds to the uncertainty 
of the initial epoch, while the latter to the offset of the orbital 
period. 

All the data were analyzed with the 6.0.0 version of the 
XMM Science Analysis System. The pipeline processing tasks 
EMCHAIN, EPCHAIN and RGSPROC were executed using the 
available calibration files, and data were filtered via EVSELECT 
to select good event patterns. We use all valid event patterns 
(PATTERN 0-12) for the MOS cameras but use valid PN events 
with single and double events (PATTERN 0-4). In the energy 
band below 1.8 keV, the emission-line peaks of the EPIC spec- 
trum are not consistent with those of the high resolution RGS 
spectrum. We adopt the RGS spectrum and ignored the EPIC 
spectrum in this energy band to avoid the statistical dominance 
by EPIC data. 

3. Results 

Figure[T]shows the X-ray image taken with MOS 1. 9 Muscae is 
the dominant X-ray source in the MOS field of view. We derived 
light curves in the soft band (0.3-2 keV) and hard band (2-10 
keV). The average count rates are 0.13 and 3.1 x 10~ 2 counts s" 1 
respectively. We detected no significant variability. 

We obtained X-ray spectra from the RGS and EPIC (PN, 
MOS 1 and 2) data. The broad band spectra are shown in 
Figure [2] and the enlarged spectrum for the RGS data is shown 
on a linear scale in Figure [3] Emission lines from helium-like 
and hydrogen-like ions of various elements (e.g. carbon, oxy- 
gen and neon) were detected. Since a single temperature plasma 
cannot reproduce these emission lines simultaneously, we ex- 
ecuted simultaneous fitting of the three datasets (RGS, MOS 
and PN) with a multi-temperature plasma model with an absorp- 
tion model (wabs) using XSPEC ver. 12. The multi-temperature 
model we adopted first is the simple collisional ionization equi- 
librium (CIE) model (cevmkl), where the emission measures fol- 
low a power-law in temperature, i.e, the emission measure from 
the electron temperature AT is proportional to {kT I kT max ) a with 
an index a ( |Singh et al.|1 996). The secondary adopted model is 
the two temperature plane-parallel shock model (2T vpshock), 
which does account for non-equilibrium ionization (NEI). 

In both fittings, line broadening and shifts were seen. Then 
we fitted the high-intensity O VHI a line around 0.65 keV with 
a single Gaussian model to investigate the width. The obtained 
broadening is O^j' 3 , eV in 1 <x, which corresponds to a velocity 
of 1 190-1720 km s 1 in FWHM. With this value, we adopted 
a Gaussian smoothing (gsmooth) of the multi-temperature mod- 
els, fixing the ratio cr/E ~ 1.3 eV / 0.65 keV and allowing the 
redshift to be a free parameter. 

The above fittings reasonably reproduced the overall spec- 
trum. The solid line in Figure [3] indicates the model for NEI. 
However, residuals around 0.49 keV were left (see Figure [3]l in 
the RGS band. The excess is identified with the RRC structure of 
C VI. We then added a recombination edge (redge) model to the 
above models, and re-fitted the datasets. The results are shown 
in Table Q] The better fit is obtained with the NEI model. The 
solid lines in Figure|2]show the best-fit NEI model. 

The best-fit NEI model, and also the CIE model, showed that 
carbon is overabundant by at least an order of magnitude com- 
pared to nitrogen in 9 Muscae (see Table [TJ, indicating that the 



plasma originates from Wolf-Rayet (WC) stellar winds. The best 
fit model has an absorption-corrected X-ray luminosity (Lx) of 
3.1 x 10 33 erg s" 1 (0 3-8.0 keV band) and log L x /L bol of -5.56 
(log L\, \/Lq = 5.47; |Nugis & Lammers|2000| , which is slightly 
larger than the known correlation Lx - 10 Lboi for O+O and 
WN+O binaries ( |Berghoefer eTal^[lW7) |Oskinova||200"5] ). The 
absorption of 2.4 x 10 zl cm~ z is consistent, within a factor of 
1.5, to the interstellar absorption 1.6 x 10 21 cm" 2 derived from 
Ay = 0.93 mag and the correlation in Predehl & Schmitt ( 1995 ). 



The red-shift is 2.2 x 10~ 3 , which corresponds to 650 km s~\ 
Independently, we further fitted each emission line with a single 
Gaussian model, and obtained a similar red-shift ~ 2 x 10~ 3 (see 
Table[2). 

The RRC structure of carbon we detected is the second ex- 
ample among WR massive binaries, with the fi rst case being 
the spectrum of y 2 Velorum ( Schild et al. 2004 1. The fitting of 
the RRC structure in 9 Muscae needed a plasma temperature of 
4.7 (3.3-6.2) eV, which is similar to but about 1.3 times higher 
than that of y 2 Velorum. The cool plasma component which 
makes up the RRC structures may be common in WR binaries. 
Interestingly, the best-fit value of the RRC indicated a red-shift 
of 2.5xl0~ 3 , which is similar to that obtained in lines, although 
the error is large. 

We further analyzed the emission lines around 0.56 keV and 
0.9 keV, which are the helium-like triplets of oxygen and neon, 
respectively. The intensity ratio of forbidden lines (HE6) to res- 
onance lines (HE4) indicates that the lines originate f rom a thin- 
thermal plasma whose density is less than 10 n cm" 3 (Porquet et 
|al.|2001) , which favors the NEI model. 

4. Discussion 



The X-ra y spectra fro m the massive bin aries WR140 ( Poll ock et 
| al.|2005| ), y 2 Velorum (jSchild et al.|2004[> WR25 ( |Raassen et al 
|2003[ >, V444 Cygnus (Maeda et al.|1999) are all strong in hard 
X-rays with high temperature components of several keV. The 
hot component of 9 Muscae also strongly suggests that X-rays 
originate from wind-wind collisions. The absorption-corrected 
X-ray luminosity of 9 Muscae is 3.1 x 10 33 erg s _1 , which is 
similar to those of the other massive binaries (WR140, 2.0 x 
10 34 erg s _1 ; y 1 Velorum, 8.4 x 10 32 erg s _1 ; WR25, 1.3 x 10 34 
erg s -1 ; V444 Cygnus, 1.4 x 10 33 erg s -1 ). 

High dispersion sp ectra of the massiv e binaries ha ve been 
obtained f or WR140 ( (Pollock et al.|2005j l, y 2 Velorum ( |Henley 
|et al. 2005 1, and 9 Muscae (this study). Significant Doppler shifts 
are only reported by |Pollock et"aL] ([2005) for WR140 at the 
pre-periastron phase when one side of the bow-shock zone ap- 
proaches toward the line of sight. They detected a significant 
blue shift of 600 km s _1 . No detection of the Doppler shift 
at post-periastron was interpreted to result from the bow-shock 
zone being aligned perpendicular to the line of sight. No detec- 
tion from y 2 Velorum is also attributed to a similar situation or 
to a wide open angle of the bow-shock zone. 

Our detection of the Doppler-shift in the emission lines is the 
second example among the massive binaries. The red-shift in the 
emission lines strongly suggests that one side of the bow-shock 
zone is receding along the line of sight, indicating that the pri- 
mary WR star is in front. However, this geometry is inconsistent 
with the orbital solution for the short -period binary in Hill et al. 
(2002) in which the O star is likely in front. Instead, if the O 
supergiant separated by 46 mas from the short-period binary is 
a companion of this system and located behind the primary WR 
star, as shown in Figure [4] the wind-wind collision zone could 
be receding. 
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Fig. 3. Unfolded XMM-Newton RGS spectrum of 9 Muscae made with S AS task rgsfluxer. The solid lines indicate wabs * gsmooth 
* 2T vpshock model. The RRC model (redge model) is not included in the model (see text). 



This idea that the wind-wind collision zone is between the 
short-period binary and the O supergiant is the same as that of 
Dougherty & Williams (2000 1. By assuming that the distance to 
6 Muscae is 2.27 kpc, the projection of 46 mas corresponds to 
~100 AU. It is well known that X-rays from massive binaries 



suffer absorption from the WR wind if the WR star is in front 
(e.g., |Maeda et al.|19 99). An X-ray emitting region widely ex- 
tended over 100 AU should avoid the absorption in the 6 Muscae 
system. This interpretation strongly suggests that the astromet- 
ric O supergiant is a companion to the system. The wind from 
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Table 1. Results of spectral fitting with the multi-temperature models. Errors and upper limits are at 90% confidence level. 



Model 



wabs * gsmooth * (redge + cevmkl) wabs * gsmooth * (redge + 2T vpshock) 



Parameter 



Best-fit values 



Best-fit values 



Line of sight absorption (wabs) 
N R (10 21 cm" 2 ) 

Doppler broadening (gsmooth) 
o- S [](eV) 

Index[] 

Doppler shift (redshift) 

Redshift (x 1(T 3 ) 

RRC (redge) 

Threshold energy (eV) 

kT RRC (eV) 
norm^x 1CT 4 ) 
Thin thermal plasma model 
a 

kT max (keV) 
kT, (keV) 
kT 2 (keV) 
Abundance^] 

C 

N 

O 

Nc 

Si 

S 

Ar 
Ca 
Fc 

t„(x 10" scirr 3 ) 
normi0(xlO" 4 ) 
norma ' (X 1(T 4 ) 



L x [|(erg s- 1 ) 



*7d.o.f 



2 56 +007 

""-0.06 



n i q+0.08 
2 - i8 -0.26 



4±1 
1.5+0.3 
(cevmkl) 
-0.25 ± 0.02 

q q+0.7 
y - y -0.5 



14.8! 34 

< io :g 

1.9+0.1 

1.8+0.3 

1.0+0.1 

2.2+0.2 
4 5+0.8 

3.7+1.1 
0.21+0.03 

4.4+0.1 

3.6 x 10 33 



11.95 (fixed) 
1 (fixed) 



o 4+0.2 



9 i 7+O.I6 
*" L ' -0.06 



489.2«;° 



1.1+0.2 
(2T vpshock) 



0.58!- 3 
3.6+0.1 



3 5 +LS 

J -0.6 

<0.1 
0.33+0.09 
47 +010 

q 3 .+0.05 

64 +010 

"•"^-0.07 

1.1+0.2 

1 o +0 - 4 

15 +002 

u - 1J -0.03 

4.9 

5.8!- 
8.3!°J 
3.1 x 10 33 



762/429 



605/426 



8 (T6 is the standard deviation at 6 keV.The broadening corresponds to a velocity of 1410 km s _I in FWHM. 
b This is the power of energy for sigma variation. 

c Normalization constant of the r edge model defined as tota l photons crrT 2 s" 1 in the line. 
d The solar abundances are from Anders & Grevesse ( 



1 1989} . For the elements with no emission lines (hydrogen, helium, sodium and aluminum), 



the solar abundances are kept fixed. The abundance of magnesium and nickel are linked to that of silicon and iron respectively. 

' Normalization constant defined as EM x 10 _14 (4^D 2 ) -1 , where EM is the emission measure in cirT 3 and D is the distance in cm. In the 
cevmkl model, emission measure follows a power law in temperature, i.e. the emission measure from electron temperature kT is proportional to 

(kT/kT max f. 

f The absorption-corrected luminosity (0.3-8.0 keV) was calculated assuming a distance of 2.27 kpc. 



the spectroscopic O star companion should be much weaker than 
that from the optically identified O supergiant. 

Another interesting possibility was obtained from the RRC 
structure. Table [2] shows that the RRC may be red-shifted with 
a similar velocity to the emission lines. If this is true, the RRC 
may originate from ions that escape from the bow-shock layer 
through diffusion or convection processes, and exchange elec- 
trons in the following less-ionized winds. Follow-up observa- 
tions with much deeper exposure are necessary to test this idea 
since we need to restrict the Doppler shift of the RRC structure 
more tightly. 



5. Summary 

The XMM-Newton spectra of 6 Muscae show He- and H-like 
emission lines from various elements as well as the RRCs from 
carbon. The results of our study are as follows: 

A. The He-like and H-like emission can be reproduced by a 
multi-temperature plasma model. A better fit was obtained 
for the NEI case (vpshock) rather than the collisional equi- 
librium (cevmkl) model. The high temperature component 
of the NEI model is as high as kT ~3 keV. The high 
temperature component indicates the existence of plasma 
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Table 2. Doppler shift of emission lines and RRC. Errors are at 90% confidence level. 



Ion" 

£„ bs (keV| 
Redshift[](x 1(T 3 ) 
v (km s -1 ) 



C VI 


C VI 


o vn 


o vm 


NelX 


NeX 


CRRC 


0.3675 


0.4356 


0.5739 


0.6536 


0.9050 


1.0218 


0.4900 


0.3668 


0.4344 


0.5725 


0.6522 


0.9022 


1.0196 


0.4892 


1 q+o.9 
1 - y -l,9 

560 +270 


2.8±0.6 


2.5±0.5 
750 +140 


2.1±0.3 

630 +8 ° 


3 l +u 

Q 30 +320 


9 9+0.8 
^•^-0.9 

650 +240 


1 6 +L8 

4gn+550 
^"-611 



Doppler shifts of the emission lines were measured with a single Gaussian fitting. That of the RRC was picked from the best-fit parameter of 
the 2T vpshock model in Table [T] 
* Theoretical energy of the line center. 
c The center energy of the observed line. 
d The redshift defined as (E - E obs )/E obs . 
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heating by a wind-wind collision shock. 

B. We detected RRC structure which implies the existence of a 
cooler component of around kT ~ 5 eV. 

C. The abundance of carbon is at least one order of magnitude 
higher than that of nitrogen. The over-abundance indicates 
that the Wolf-Rayet (WC) stellar winds dominate the X-ray 
emitting gas. 
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D. The emission lines from carbon, oxygen, neon and silicon 
and possibly the RRC from carbon show red-shifts of ~ 600 
km s~' with broadenings of ~ 1400 km s -1 in FWHM. The 
red-shift would be evidence supporting the widely separated 
O supergiant as a companion with which the collision zone 
could be formed, lying behind the short-period binary. The 
wind from the O star in the short-period binary should be 
much weaker than that from the O supergiant in the wide 
binary. 



